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We report the results ol theoretical investigations of low frequency tunneling current noise spectra 
component (l// a ). Localized states of individual impurity atoms play the key role in low frequency 
tunneling current noise formation. It is found that switching "on" and "off' of Coulomb interaction 
of conduction electrons with one or two charged localized states results in power law singularity of 
low-frequency tunneling current noise spectrum 1/ f a . Power law exponent in different low frequency 
ranges depends on the relative values of Coulomb interaction of conduction electrons with different 
charged impurities. 
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I. INTRODUCTION 

In the present work we discuss one of the possible rea- 
sons for the l/f a noise in the STM/STS junctions. We 
suggest the theoretical model for tunneling current noise 
above the flat surface and above the impurity atoms on 
semiconductor or metallic surfaces. We found out that 
changing of the power law exponent above the flat surface 
and above the impurity atom depends on the parameters 
of tunneling junction such as tip-sample separation 

Problem of low frequency noise with 1//" spectra for- 
mation in electron devices is one of the most interest- 
ing and important in recent years. Usually the typical 
approach to l/f a noise problem consists of "by hand" 
introducing of the random relaxation time To for two- 
state system with the probability distribution function 
A/tq. Therefore the averaged over tq noise spectra of 
two-states system has power law singularity. But the 
physical nature and the microscopic origin of such prob- 
ability distribution function in general is unknown. Up to 
now only the limited number of works was devoted to 
the problem of 1//" noise study. The investigations of 
the noise in two-level system was carried out in 1. Au- 
thors studied current noise in a double-barrier resonant- 
tunneling structure due to dynamic defects that switch 
states because of their interaction with a thermal bath. 
The time fluctuations of the resonant level result in low- 
frequency noise, the characteristics of which depend on 
the relative strengths of the electron escape rate and the 
defect's switching rate. If the number of defects is large, 
the noise is of the 1// type. In [j| authors studied shot 
noise in a mesoscopic quantum resistor. They found cor- 
relation functions of all order, distribution function of 
the transmitted charge and considered Pauli principle 
as the reason for the fluctuations. Altshuler et al. Q 
studied current fluctuations in a mesoscopic conductor. 
They derived a general expression for the fluctuations 
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in the cylindrical tunneling contact in the presence of 
a time dependent voltage. In [5[ Moller, Esslinger and 
Koslowski have investigated the noise of tunneling cur- 
rent at zero bias voltage. The measurements were carried 
out at UHV conditions at base pressure 5 x 1CP 11 torr. 
The authors have demonstrated that at zero bias voltage 
the 1/ f a component of noise in the tunneling current 
vanishes and white noise becomes dominant. Tiedje et 
al. @ have found the 1//" dependence of the current 
noise in STM experiments on graphite in ambient condi- 
tions. They attributed this effect to fluctuations induced 
by adsorbates in tunneling junction area. 

In [?J the fluctuations of tunneling barrier height have 
been investigated. The experiments have been performed 
under UHV conditions on graphite and gold samples us- 
ing Ptlr tips. From these measurements the authors have 
concluded that the intensity of barrier height fluctuations 
correspond to the intensity of tunneling current l/f a 
noise in the frequency range from 1 to 100 Hz. 

Our aim is to study one of the possible microscopic ori- 
gins of 1// Q noise in tunneling contact. We shall anal- 
yse tunneling contact model with localized states and 
Coulomb interaction between them which gives us op- 
portunity to describe the formation of 1//" tunneling 
current noise spectra. 



II. THE SUGGESTED MODEL AND MAIN 
RESULTS 

Let's start from the model of two localized states in 
tunneling contact. In this case one of the localized states 
is formed by impurity atom in semiconductor and the 
other one by tip apex. 

When electron tunnels in or from localized state, the 
electron filling numbers of localized state rapidly change 
leading to appearance of localized state additional charge 
and sudden switching " on" and " off' Coulomb potential. 
Electrons in the leads feel this Coulomb potential. 

The model system (Fig. 1) can be described by hamil- 
tonian H: 



2 




E F -eV 



FIG. 1: Schematic diagram of tunneling processes through 
states localized on impurity atom and on the STM tip apex. 
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Hq describes free electrons in the leads and in localized 
states. H tU n describes tunneling transitions between the 
leads through localized states. H int corresponds to pro- 
cesses of intraband scattering caused by Coulomb poten- 
tials Wi, Wi of localized states charges. 

Operators c k (c k ) and c+ (c p ) correspond to electrons in 
the leads and operators af(a,i) correspond to electrons in 
the localized states with energy e,. 

Current noise correlation function is determined as: 



(H/e) 2 ■ S(t,t r ) =<I L (t)-I L (t r ) >= 

J2 n<c+{t')a t {t')a+{t)ci{t)> 
k,k 



(2) 



where 



$>* • e = 4(t') ai (t')T ki - h.c.) ■ - (3) 

k k 



The current noise spectra is determined by Fourier trans- 
formation of S(t,t'): S(uj) = J S{r)dT ■ e luJT . 

We shall use Keldysh diagram technique in our study 
of low frequency tunneling current noise spectra [8[. 
Tunneling current noise spectra S(u>) can be expressed 
through Keldysh Green functions. 

For the one localized state in tunneling contact one 
should put: T = 0, T k2 = 0, e 2 = 0. 

Expression for tunneling current noise spectra without 
Coulomb re-normalization of the tunneling vertexes can 
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FIG. 2: Lowest order diagrams contributing to tunneling cur- 
rent noise spectra for one localized state in tunneling contact, 
a) In the absence of Coulomb re-normalization of tunneling 
vertexes. b) In the presence of Coulomb re-normalization of 
tunneling vertexes - ladder diagrams, c) In the presence of 
Coulomb re-normalization of tunneling vertexes - parquet di- 
agrams. Tunneling vertexes are marked by the dot. Coulomb 
potential is marked by the bold dot. Solid lines correspond 
to electron Green functions. 



be found from diagrams shown on Fig. 2a and has the 
form: 



h /ef • S (w) = T 2 kl • ]T(G< fci H ■ G> [u + J) + 



k.k± 

+ G<H • G> ki (w + a/)) + T 2 kl ■ £(G<( 

■ G> 1 (c + c') + G< 1 M-G>( w + ^)) (4) 
Green functions are evaluated from equations: 

^fefc! - W'fc! + ^kk 2_, I k' U k ' kl + U kk 2^1 1 k' hr k 'k 1 



kl — U kk J fel ( - r ll + ^fefe J-kl^H 



(5) 



G< 



-lin (w)/mGf 1 (w) 



Functions G kik and G{* fe can be found by substitution 
ni on iii ^ 1- 



ni(u>) 



7fci + 7pi 

are non-equilibrium localized state filling numbers. 



(6) 
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n£(w) and n°(co)- equilibrium filling numbers in the 
leads. In our model relaxation rates 7^, 7^ are deter- 
mined by electron tunneling transitions from localized 
states to the leads k and p continuum states: 
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rp2 /^/-(O-R _ . 
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TllmGl* = 7fe 



(7) 



For tunneling current noise spectra without Coulomb re- 
normalization of tunneling vertexes after substitution the 
correspondent Green functions can be written as: 

( h/e) 2 ■ S (uj) = Tfei ■ / dw'-TmGfiCu') • 

• ImGf x {u + u') ■ (m(u} + u}') - 1) • 

• (ni(w') - n fc (a/)) + ni(a/) • {n\{uj + ui') - 

— rife (a; + u') + 7^ • / rfcj'7mG r f 1 (a;') • 

• ImGii(u + J) ■ (nk{u + u') - 1) • ni(w') - 

— ni(w') • (ni(w + U)') — 1) — n fe (aj + a/) • n k {uj') + 

+ rtfe(w') • (tt-! ((jj + cl>') + 7fc! • J dui 1 'ImGf 1 (w + a/) • 

• (n fc (t/)) • (m(w + w') - 1) + JroGgV) ■ 

• (n 1 (w / ))-(n fc (w + w , )-l)=5o (8) 

This expression gives us an opportunity to analyze tun- 
neling current noise spectra for typical values of kinetic 
parameters of tunneling contact when localized charge is 
connected with the tip apex state eV — £\. Some low 
frequency spectra are shown on Fig. 3a. 

It is clearly evident that when frequency aspire to zero 
tunneling current spectra aspire to constant value. 

Now let us consider re-normalization of the tunneling 
amplitude and vertex corrections to the tunneling current 
noise spectra caused by Coulomb interaction between 
charged localized state and the electrons in tunneling 
contact leads. The result of re-normalization is shown on 
Fig. 2b and Fig. 2c. Re-normalization gives us two types 
of diagrams contributing to the final tunneling current 
noise spectra expression. Ladder diagrams is the most 
simple type of diagrams which gives logarithmic correc- 
tions to vertexes. Ladder diagrams give logarithmic di- 
vergency at the threshold voltage eV — s\ (Fig. 2b). 
But this is not the only relevant kind of graphs. We 
must consider one more type of graphs (parquet graphs) 
which gives logarithmically large contribution to tunnel- 
ing spectra when u — > and eV = E\ (Fig. 2c). In par- 
quet graphs a new type of "bubble" appears instead of 
"dots" in ladder diagrams. In this situation one should 
retain in the n-th order of perturbation expansion the 
most divergent terms. For the first time this method was 
developed by Dyatlov et. al. [l(|. It was shown that 
for proper treatment of this problem one should write 
down integral equations for the so called parquet graphs, 
which are constructed by successive substitution the sim- 
ple Coulomb vertex for the two types of bubbles in per- 
turbation series. The integral equations can be solved 



with logarithmic accuracy, as it was done, for example, 
by Nozieres for edge singularities in X-ray absorption 
spectra in metals 11| . The expression for tunneling cur- 
rent noise spectra after Coulomb re-normalization can be 
written as: 



( H/e) 2 -S(cj) = Sv(uj) + S {lo) 
D 2 



(( 



(LU + eV-ex) 2 + (7fci +7 P i) 2 ' 

D 2 

■ (-w + eV - £i) 2 + (7fci +7 P i) 2 



) Wv ) (9) 



where D is the bandwidth for electrons in tunneling con- 
tact leads, , W-Coulomb potential, v- the equilibrium 
density of states in the tunneling contact leads. We con- 
sider that localized state is formed by STM tip apex, so 
we have to put eV = E\, Then we obtain final expression: 



(fi/e) 2 -5H = S*oH(l + (- 



D 2 



) Wv ) (10) 



+ ilki +l P if 



Fig. 3b-3c demonstrate low frequency tunneling cur- 
rent noise spectra for typical values of dimensionless ki- 
netic parameters. We can see that re-normalization of 
tunneling matrix element by switched "on" and "off" 
Coulomb potential of charged impurity lead to typical 
power law singularity in low frequency tunneling current 
noise spectra. 

Some tunneling current noise spectra in double loga- 
rithmic scale for typical values of kinetic parameters are 
shown on Fig. 3d-3e. Low frequency tunneling current 
noise spectra in double logarithmic scale make it clear 
that power law exponent depends only on Coulomb po- 
tential of charged impurity and does not depend on tun- 
neling contact parameters. 

Now let's describe interaction effects of two local- 
ized states in tunneling contact: surface localized state, 
formed by impurity atom, and localized STM tip apex 
state (Fig. 1). Expression which describes tunneling cur- 
rent noise spectra without Coulomb re-normalization can 
be calculated from graphs, shown in Fig. 4a. It consists 
of three parts. 



Sq{oj) = S 01 (w) + #02 M + S 03 (u}) 



(11) 



Sai(u>) and So2(w) are rather simple parts equal to 
tunneling current noise spectra without Coulomb re- 
normalization in the casejff one localized state in tun- 
neling contact Eq. ([5]). S03 ( w ) is not trivial, it exists 
only due to electron tunneling transitions from one lead 
to both localized states. Green functions shown on the 
graphs are found in The contribution of S'oi(w) is 

given by graphs with 1 = 7 = 1 (Fig. 4a), £02 (w) is de- 
scribed by graph with % — j = 2 (Fig. 4a), So3(w) is given 
by diagrams with i ^ j (Fig. 4a). The final expression for 
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FIG. 3: Typical low frequency tunneling current noise spectra for different values of dimensionless kinetic parameters for 
one localized state in tunneling contact (eV = ei = 1). a) In the absence of Coulomb re-normalization of tunneling vertexes. 
b)-c) In the presence of Coulomb re-normalization of tunneling vertexes. d)-e) In the presence of Coulomb re-normalization of 
tunneling vertexes in double logarithmic scale. 



tunneling current noise spectra for two localized states in 
tunneling contact without Coulomb re-normalization of 
tunneling vertexes after substitution the correspondent 
Green functions can be written as: 



( hie) 2 ■ S 03 (cj) = 8 • 7fcl • lk2 ■ J du'ImGf^u') • 

■ ImG 22 {uj + u>') ■ (ni(w') • (n 2 (w + o/) - 1) + 

+ n k (ui') ■ (n 2 (ui + ui') - 1) + (m(u/) • 

• (n k (uj + w 1 ) - 1) - n k (uj') ■ (n k (uj + w 1 ) - 1) + 

+ 8 • 7fcl • 7fe2 ■ / du'ImG^uo 1 ) • JmGf^w + w') • 



• (n 2 (u') ■ + u)') - 1) + 

+ n k (oj') ■ (ni(w + oj') - 1) + {n 2 {uj') ■ 
■ (n*(w + <J) - 1) - n k {u') • (n*(w + w') - 1) (12) 



Non-equilibrium filling numbers ni and n2 are deter- 
mined from Dyson equations for Keldysh functions G*j-, 
where i,j = 1,2 in [12]: 

Green functions G{\ , G 22 have rather simple form and 
have been also derived in 12]. We don't take in ac- 
count diagrams including Green functions G^ 2 ,G 21 . In 
our case of weak interaction between localized states 
(T < 7fci,7fe2,7p2) contribution to the noise spectra of 
such diagrams has additional small parameter in com- 
parison with diagrams depicted on Fig. 4a: 



T ■ 7fei 



(ei - e 2 ) 2 + (7fei + 7fc2 + 7 P 2) 2 



T/ 7 



(13) 



Some typical low frequency tunneling current noise spec- 
tra for different values of dimensionless kinetic param- 
eters without Coulomb re normalization are shown on 
Fig. 5a. It is clearly evident that when frequency aspire 
to zero tunneling current spectra aspire to constant value 
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FIG. 4: Lowest order diagrams contributing to tunneling cur- 
rent noise spectra for two localized states in tunneling contact, 
a) In the absence of Coulomb re-normalization of tunneling 
vertexes. b) In the presence of Coulomb re-normalization of 
tunneling vertexes. Tunneling vertexes are marked by the dot. 
Coulomb potential W\ is marked by the bold dot. Coulomb 
potential W2 is marked by the star. Solid lines correspond to 
electron Green functions. 



for different dimcnsionlcss kinetic parameters. Now let us 
consider re-normalization of the tunneling amplitude and 
vertex corrections to the tunneling current spectra caused 
by Coulomb interaction between both localized states 
and tunneling contact leads. Re-normalization gives us 
two types of diagrams contributing to the final tunneling 
current noise spectra expression similar to types of graphs 
in the case of one localized state in tunneling contact. It 
is necessary to re-normalize each vertex individually and 
re-normalize both vertexes jointly. For each diagram re- 
normalized by Coulomb potential of localized state one 
have to take in to account the whole series of diagrams in- 
cluding Coulomb interaction with another localized state 
(Fig. 4b). 

The final expression for tunneling current noise spectra 
after Coulomb re-normalization of tunneling vertexes can 
be written as: 



( h/e) 2 -S(cj) = S (u) + S 01 (cj)- 

D 2 , D 2 



(( 



+ S02(w)-(( 



D 2 



{u + Etf + T 2 ' \u J + E 2 ) 2 + T 2 l 



+ ( 



D 2 



) w n + s Q3 (u) ■ 



(14) 



D 2 



D 2 
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{u + E 1 y + T 2 ' > + £ 2 ) 2 + r 2 

D 2 , w .„ , D 2 



-{-uj + E^ + T 2 ' ' '(-cu + E 2 ) 2 +r 
where, 



£1 



zV - 



si+e 2 -e 2 ) 2 +4T 2 



E 2 = eV-^-^- E f +4T ' 2 (15) 



r 2 



2 2 

(gi - £2)(7fc2 + 7p2 - 7fci) , 7fc2 + 7p2 - Tfci 
2V(£i-e 2 ) 2 +4T 2 r 2 

(£1 - £ 2 )(7fc2 + 7 P 2 - Tfci) _ 7fc2 + 7 P 2 - 7fci 
2 V / (£i-£ 2 ) 2 +4T 2 ~ 2 



When the applied bias voltage become close to e 2 we 
obtain low frequency tunneling spectra. So we have to 
put eV = £1 = e 2 . We obtain the expression: 



( h/e) 2 ■ sh = SoH + S01H • {{-^r^ + 



+ ( 



D 2 , 
D 2 



+ 5o3(w)-((- 



•^ 2 + r 2 ' 

D 2 



) WlV + (16) 



LO 



2 + r 2 



LO 



2 + Tf 



LO 



2 + r 2 - 



; 2 + r 2 - 



First of all we consider the situation when both lo- 
calized states acquire positive charge. Fig. 5b demon- 
strate low frequency tunneling current noise spectra for 
typical values of dimcnsionlcss kinetic parameters. We 
can see that re-normalization of tunneling matrix ele- 
ment by switched "on" and "off" Coulomb potential of 
charged impurities lead to typical power law singularity 
in low frequency tunneling current noise spectra. Tun- 
neling current noise spectra in double logarithmic scale 
demonstrate frequency regions where every part of final 
expression which include different power law exponent, 
approximate noise spectra in the best way (Fig. 5c). Let's 
analyse tunneling current spectra shown on Fig. 5c. Con- 
tribution to tunneling current spectra at low (zero) fre- 
quency is always determined by the term which produces 
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FIG. 5: Typical low frequency tunneling current noise spec- 
tra for different values of dimensionless kinetic parameters for 
two localized states in tunneling contact (eV = ei = £2 = 1). 
a) In the absence of Coulomb re-normalization of tunneling 
vertexes. b) In the presence of Coulomb re-normalization 
of tunneling vertexes. c) In the presence of Coulomb re- 
normalization of tunneling vertexes in double logarithmic 
scale. Black lines correspond to tunneling current noise spec- 
tra in double logarithmic scale. Gray lines correspond to the 
terms which produce the most strong logarithmic singularity 
and approximate full spectra in the best way. 



the most strong of logarithmic singularity, determined by 
the sum of localized states Coulomb potentials at any 
values of dimensionless tunneling rates of tunneling con- 
tact. In the case of W\ » W2 (Coulomb potential of 
charged impurity atom is larger than Coulomb potential 
of charged tip apex localized state) with frequency in- 
creasing the main contribution to tunneling current noise 
spectra is given by the term depending on Coulomb po- 
tential of localized state caused by impurity atom (Wi) 
at any dimensionless kinetic parameters of tunneling con- 
tact. If Coulomb potential of charged tip apex localized 
state is larger than Coulomb potential of charged impu- 
rity atom (W2 >> W\) with frequency increasing the 
main contribution to tunneling current noise spectra is 



given by the term depending on Coulomb potential of 
tip apex localized state (W2) at any dimensionless ki- 
netic parameters of tunneling contact. When Coulomb 
potential of charged tip apex localized state is similar to 
Coulomb potential of charged impurity atom (W\ ~ W2) 
with frequency increasing the main contribution to tun- 
neling current noise spectra can be given both by the 
term depending on Coulomb potential of tip apex local- 
ized state (W2) and by the item depending on Coulomb 
potential of localized state caused by impurity atom (Wi) 
at any dimensionless kinetic parameters of tunneling con- 
tact. In this case the replacement of dominating term 
can take place. Simple estimation gives an opportunity 
to determine the validity of obtained expressions. For 
typical typical parameters of tunneling junction model 
in the zero frequency region power spectrum of tunnel- 
ing current corresponds to experimental results. Power 
spectrum on zero frequency has the form: 



5(0) « ( 7 e// ■ e/hf ■ {D/ lein r w ■ (1/ A w ) (17) 
5(0) « 10~ 18 A 2 '/Hz [J. 



For typical y efJ , 7e m « 10~ 13 , D « 10, W w 0, 5 



III. CONCLUSION 

The microscopic theoretical approach describing tun- 
neling current noise spectra taking in account many- 
particle interaction was proposed. When electron tun- 
nels to or from localized state the charge of localized 
state rapidly changes. This results in sudden switching 
on and off of additional Coulomb potential in tunnel- 
ing junction area, and leads to typical power law depen- 
dence for low frequency tunneling current noise spectra. 
In the case of two interacting positively charged local- 
ized states the power law exponent at low frequency is 
determined by the sum of Coulomb potentials. It was 
shown that if one of the Coulomb potentials strongly ex- 
ceeds the other one this potential determine tunneling 
current noise spectra at low frequency with increasing of 
frequency. If W\ ~ W2, tunneling current noise spec- 
tra can be determined both by Coulomb potential Wi, 
and by Coulomb potential of charged tip apex state W2 , 
depending on dimensionless kinetic parameters of tun- 
neling contact. When localized states acquire charges of 
opposite signs tunneling current noise spectra in the low 
frequency region are approximated by the term depend- 
ing on maximum positive value of Coulomb potentials 
Wi, W 2 , Wi - W 2 or W 2 - Wi. 
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